arXiv: 1505.00610v2 [math.PR] 3 Sep 2015 


Correlation kernels for sums and products of 
random matrices 
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Abstract 

Let A be a random matrix whose squared singular value density 
is a polynomial ensemble. We derive double contour integral formulas 
for the correlation kernels of the squared singular values of GX and 
TA, where G is a complex Ginibre matrix and T is a truncated uni¬ 
tary matrix. We also consider the product of A and several complex 
Ginibre/truncated unitary matrices. As an application, we derive the 
precise condition for the squared singular values of the product of sev¬ 
eral truncated unitary matrices to follow a polynomial ensemble. We 
also consider the sum H + M where iJ is a GUE matrix and M is 
a random matrix whose eigenvalue density is a polynomial ensemble. 
We show that the eigenvalues oi H + M follow a polynomial ensemble 
whose correlation kernel can be expressed as a double contour integral. 
As an application, we point out a connection to the two-matrix model. 


1 Introduction 

Eigenvalues and singular values of sums and products of random matrices 
have been studied extensively over the last decade. The focus of the research 
has been the global density of the eigenvalues and singular values of the sum 
and product of large random matrices. As the dimension of random matrices 
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tends to infinity, the macroscopic eigenvalue density of sums and products 
of random matrices can be described under fairly general conditions using 
free probability techniques m- However, the free probability method can 
only give global asymptotic results. For exact correlation functions of the 
eigenvalues and singular values of the sums and products of random matri¬ 
ces, until recently there were only a limited number of cases in which results 
were known. Notable cases are the Gaussian unitary ensemble with exter¬ 
nal source [9], which is the sum of a Gaussian unitary matrix (GUE) with 
a deterministic matrix; and the complex Wishart ensemble [S], which is the 
product of a complex Ginibre matrix with a deterministic matrix. Another 
interesting result is known for the sum of complex null Wishart matrices, in 
the form of the multiple Laguerre minor process with a fixed time [T]. It is 
remarkable that all the examples mentioned above are polynomial ensembles, 
a special kind of determinantal point processes, and that their correlation 
kernels all have double contour integral representations. 

Recently, the correlation functions of the (squared) singular values are 
found for more types of products of random matrices, for instance for prod¬ 
ucts of complex Ginibre matrices and products of truncated unitary Haar 
distributed matrices, as well as their inverses, see [3], [28], HZ], [27], m, 
and [121 . For a comprehensive survey of the current developments, see [2] 
and references therein. In all the cases listed above, the squared singular 
value densities are polynomial ensembles, and the correlation kernels can be 
expressed as double contour integrals. 

Polynomial ensembles arise naturally in the study of probability densi¬ 
ties for eigenvalues and singular values of random matrices. Their structure 
is preserved under certain operations on random matrices, such as multipli¬ 
cation by a complex Ginibre matrix m and multiplication by a truncated 
unitary Haar distributed matrix [25] . Those random matrix operations thus 
induce transformations of polynomial ensembles [2B|. In this paper, we ob¬ 
tain explicit transformation formulas for the correlation kernels for a number 
of such transformations. 

One feature of our results is that the correlation kernels obtained in this 
paper are all in the double contour integral form. Double contour integral 
formulas are very suitable to derive asymptotic results on local statistics, by 
the classical saddle point method. Numerous applications of double contour 
formulas can be found in literature. Two recent examples are |28j for the 
discovery of a family of new hard-edge universality and [29] for the proof of 
bulk and soft-edge universality. 

A polynomial ensemble m is a probability density function for n parti¬ 
cles xi,... ,Xn on the real line of the form 

^A„(x)det[/fc_i(xj)]”^i (1.1) 
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with given functions /o,..., fn-i- Here we use 


^n{x) = n - Xj) = det ^ 


n 


j,k=l 


( 1 . 2 ) 




to denote the Vandermonde determinant for the n-tuple x = (xi,..., x„) G 
and Zn is a normalization constant. It is a polynomial ensemble on 
E C M if fk{x) = 0 for every x £ M. \ E and every k = 0,..., n — 1, since 
then the particles are in E (with probability one). 

The polynomial ensemble (II.ip is a determinantal point process and its 
correlation kernel can be written in the form 


n—1 



(1.3) 


j=0 


where pj is a polynomial of degree j and qj belongs to the linear span of 
/o,..., fn-i for every j = 0,..., n — 1, in such a way that they satisfy the 
biorthogonality conditions 



We call such a set of pj's and g/s a biorthogonal system associated with the 
polynomial ensemble m- We sometimes refer to the qj ’s as dual functions. 
The pj’s and qj's are not unique, but the sum Pj(^)9i(y) fo (11-311 is. 

If we require that pj is a monic polynomial, and that qj is in the linear span 
of fo,..., fj for every j, then the biorthogonal system is unique. However, 
such a biorthogonal system does not always exist, and, if it exists, it may not 
be the most convenient one to work with. The unique monic polynomial pn 
of degree n which is orthogonal to /o,..., fn-i, he. J^^Pn{x) fk{x)dx = 0 
for /c = 0,... , n — 1, is the average characteristic polynomial 



(1.5) 


where E denotes the average over particles xi,...,x,i in the polynomial 
ensemble m- 

In [2B] a number of transformations which preserve the structure of a 
polynomial ensemble are given with transformation formulas for the func¬ 
tions fk in (|1.1|) . but without explicit formulas for the correlation kernels. 
In this paper, we give transformation formulas for the correlation kernels, 
biorthogonal systems and average characteristic polynomials for three such 
transformations. The polynomial ensembles which we consider are joint 
probability densities for; 


3 




(i) the eigenvalues of the sum H+M of a GUE matrix H with a Hermitian 
random matrix M, 

(ii) the squared singular values of the product GX of a complex Ginibre 
matrix G with a random matrix X, 

(iii) the squared singular values of the product TX of a, truncation of a 
Haar distributed unitary matrix with a random matrix X. 


Here and throughout the paper, we understand M and X as follows. 

Assumption 1.1. M is an n x n Hermitian random matrix whose eigenvalues 
xi,. .. ,Xn follow a polynomial ensemble expressed by dni). We let Pkix) 
and qk{y) {k = 0,1,... ,n — 1) be a biorthogonal system, kn{x,y) be the 
correlation kernel, and Pn{x) be the average characteristic polynomial for 
the polynomial ensemble. 

Assumption 1.2. X is an £ x n {i > n) complex matrix whose squared 
singular values xi,... ,Xn, that is, eigenvalues of X*X, follow a polynomial 
ensemble expressed by We let pk{x) and qkiv) (A: = 0,1,... , n — 1) 

be the biorthogonal system, kn{x,y) be the correlation kernel, and Pn{x) be 
the average characteristic polynomial for the polynomial ensemble. 


Note that although notations fk, Pk, qk, kn and have different mean¬ 
ings in Assumptions 11.11 and 11.21 since M and X will not appear together 
in our paper, no confusion will occur. 

transformation formulas were obtained in m and 


In cases (ii) and (iii) 


for the joint probability densities of the squared singular values of GX 
and TX in terms of the density for X. We will prove a similar transformation 
formula in case |(i)l Our main focus is on the correlation kernels. We will 
the eigenvalue correlation kernel of H -\- M in terms 


express, m case 


of the eigenvalue correlation kernel of M, and in case (ii) (resp. case (iii)) 


the squared singular value correlation kernels of GX (resp. TX), in terms 
of the squared singular value correlation kernel of X. In addition, we will 
give transformation formulas for biorthogonal systems and for the average 
characteristic polynomials. 


2 Statement of main results 

2.1 Addition of a GUE matrix 

Up to an overall scale factor, a GUE matrix is a Hermitian random matrix H 
with diagonal entries in standard real normal distribution, upper-diagonal 
entries in standard complex normal distribution, and all upper-triangular 
and diagonal entries independent. The joint probability density function for 
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the eigenvalues xi,..., oi H is given by 

- n ^2 n 

—A„(x)^ = (27r)^ J|/c!, (2.1) 

]=1 k=l 

which is the polynomial ensemble dni) in the case where fk{x) = x^e 
If we add a GUE matrix to a random matrix with eigenvalues in a poly¬ 
nomial ensemble, we are led to a transformed polynomial ensemble for the 
eigenvalues of the sum. 

Theorem 2.1. Let H be an nxn GUE matrix and M be the random matrix 
from Assumption \1.1[ independent of H. The density of the eigenvalues 
2 / 1 , ■ ■ ■, yn of H + M is given by 

-^An{y) det , (2.2) 

for some constant Z!^, where 

/ oo 

fkit)e~ 2 ^y-^'> dt, k = 0,...,n-l. (2.3) 

-OO 

We will prove Theorem 12.11 in Section 13.11 

Remark 2.2. The functions in the transformed polynomial ensemble (|2.2p 
are convolutions of the functions fk from the original polynomial ensemble 
with the Gaussian density e ^ 

We now describe the effect of this transformation on the biorthogonal 
system, on the correlation kernel, and on the average characteristic poly¬ 
nomial. Since the eigenvalue density oi H M is the polynomial ensemble 
( 1221 ), the associated correlation kernel Kn can be written in the form 


n—1 

Kn{x,y) = '^Pk{x)Qk{y), (2.4) 

k=0 

where Pk is a polynomial of degree k and Qk is in the linear span of 
Fq, ..., Fn-i, such that 

/ OO 

Pjix)Qk{x)dx = 6j^k, for j. A: = 0,... ,n - 1. (2.5) 

-OO 

We express the quantities Kn, Pk, and Qk in the transformed polynomial 
ensemble in terms of their counterparts kn, Pk, and qk in the original poly¬ 
nomial ensemble for M. 

Theorem 2.3. Let H be an nxn GUE random matrix and M be the random 
matrix from Assumvtion \l. 11 independent of H, as in Theorem \2.1\ Then, 
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(a) a biorthogonal system in the transformed polynomial ensemble defined 
by (|2.2p and ()2.3p is given by 

1 p+ioo 

Pk{x) = / pk{s)e2^--^^"ds, (2.6) 

V ZTTZ J —ico 

1 

Qk{y) = -j= / qk{t)e-^^y-^'>"dt, (2.7) 

V J — oo 

(b) the correlation kernel for the eigenvalues of H + M is given by 

1 p-\-ioo poo 

Kn{x,y) = — ds dtknis,t)e-2^^^-^^"-^y-^'>"\ (2.8) 

ATTI j _2oo j — oo 

("cj the average characteristic polynomial Pn of H + M is given by 
with k = n. 

The proof of Theorem 12.81 will be given in Section [8.21 

Remark 2.4. In the integral formula (|2.8p . we use an extension of the ker¬ 
nel kn{s,t) for s in the complex plane. We can do this because kn{s,t) is 
polynomial in s for any t E M. 

Remark 2.5. The integral formulas (12.71) and (j2.8p for and Fk can be 
recognized as Weierstrass transforms [22l Chapter VIII] of qk and fk (up to 
the prefactors; the usual definition of the Weierstrass transform also has 1/4 
instead of 1/2 in the exponent); Pk given in (12.6p is the inverse Weierstrass 
transform of pk, up to the prefactor, which is chosen such that Pk is monic 
if pk is monic. 

As an application of Theorem 12.81 we derive in Section 18.81 a double 
contour integral formula (18.221) for the correlation kernel of the eigenvalues 
of the sum of a CUE matrix and an arbitrary random unitary invariant 
matrix. This formula is suitable for the analysis of limiting local statistics, 
and we investigate the analytic consequences in a separate paper. 

2.2 Multiplication with a Ginibre matrix 

A complex Ginibre matrix is a rectangular matrix whose entries are inde¬ 
pendent and identically distributed in the standard complex normal distri¬ 
bution. We define the complex Ginibre matrix G as follows. 

Assumption 2.6. Let n > 0 and i > n. G is an (n + u) x i complex Ginibre 
matrix. 

In the case i = n, the joint probability density for the squared singular 
values xi,..., of G is given by 

-^Anixfllx-ie-^fi X, >0, Zn = llk\T{k + u). (2.9) 
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This is a polynomial ensemble on [0,oo) defined by the functions fk{x) = 
/c = 0,1,..., n — 1. 

Left multiplication of a random matrix by a complex Ginibre matrix 
induces a transformation of polynomial ensembles on [0, oo): the following 
result was proved by Kuijlaars and Stivigny [27] (based on ideas from |3lll]). 

Proposition 2.7. Let G be the complex Ginibre matrix from Assumption 
1^.61 and X be the random matrix from Assumption li.^L independent of G. 
Then the squared singular values yi,... ,yn ofY = GX follow a polynomial 
ensemble on [0, oo) given by 

-^Aniy) det[Fk-i{yj)]lf.^^, ( 2 . 10 ) 

for some constant Z'^, where 

Fk{y) = Fe-^fk {^)j, y>0. (2.11) 

We complement the above result with transformation formulas for the 
biorthogonal system, for the correlation kernel, and for the average charac¬ 
teristic polynomial. 

Theorem 2.8. Let G be the complex Ginibre matrix from Assumvtion \2.6\ 
and X be the random matrix from Assumvtion M .A independent of G, as in 
Proposition |g. 7| . Then, 

(a) a biorthogonal system in the transformed polynomial ensemble defined 
by (I2.10p and (|2.11li is given by 

Qk{y)= Fe-^qk(^)j, (2.13) 

where Yi is a simple counter-clockwise oriented contour encircling the 
origin, 


(b) the correlation kernel Kn for the squared singular values 
given by 


of GX is 


Kn{x,y) 


2-kiJy,sJq t \s tJ 



(2.14) 


(c) the average characteristic polynomial Pn of X*G*GX is equal to 


Pnix) 


(n + z/)! 
27ri 


s ‘'e^'pn 



ds 

s 


(2.15) 
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The proof of Theorem 12.81 will be given in Section 15.11 

Remark 2.9. The functions and can be recognized as Mellin convo¬ 
lutions of fk and qk with the Gamma density Fe~^. We will show that 
Pfc defined in ()2.12D can also be written in an alternative way: if Pk{x) = 
Yl^=o , then we have 






(2.16) 


Thus the factor (n + z^)! in (I2.15p makes Pn{x) monic.. The transformation 
formula (I2.16p for is reminiscent of similar transforms in [181 Formulas 
(3.16) and (3.43)]. We can also write (|2.14p as 


Kn(x,y) = 


1 


27rf \x 


ds 


dt 


s •5 Jo t 


knis,t) 


X_y 

e« *. 


(2.17) 


2.3 Multiplication with a truncated unitary matrix 

The truncated unitary matrix T is defined as follows. 

Assumption 2.10. Let u>l, m>l>n and iJL = m — n — u> 1. Let U 
be an m X m Haar distributed random unitary matrix. T is the (n + n) x i 
truncation of U. 

In the case where i = n and p > n, the joint probability density for the 
squared singular values xi,..., of T is given by [T6l Section 3.8.3] 


1 ” 


\m—2n—u 


0<Xj < 1, 


i=i 




k=l 


T{k + m — n) 


The following result is due to Kieburg, Kuijlaars and Stivigny 

Proposition 2.11. LetT be the truncated unitary matrix from Assumption 
\2.1(]\ and X be the random matrix from Assumution \1.A independent of T. 
Then the squared singular values yi,... ,yn ofY = TX follow a polynomial 
ensemble on [0, oo) given by 


1 


Xn{y) det[Pfc-i(yi)]"fe=i, 


for some constant Z'^, where 

'y\ dt 


(2.19) 


niy) = f f'(i-0'-'A(f)T' 


y > 0. 


( 2 . 20 ) 














We have the counterpart of Theorem 12.81 as follows. 


Theorem 2.12. Let T be the truncated unitary matrix from Assumption 
\2.10\ a.nd X he the random matrix from A ssumvtion I i. gl independent ofT, 
as in Provosition \2. IR Then, 

(a) a biorthogonal system in the transformed polynomial ensemble defined 
by (I2.19P and (12.201) is given by 

Qk{y) = ( 2 . 22 ) 

where is a simple counter-clockwise contour around 0 but not con¬ 
taining 1, 


(b) the correlation kernel for the squared singular values ofTX is given 
by 


Kn{x,y) 



(2.23) 


(c) the average characteristic polynomial Pn of X*T*TX is equal to 


Pn{x) 


/i!(n + v)] 
2'Ki{n + u + fj,)\ 



_ s)->^-^Pn 



ds 

s 


(2.24) 


We will prove Theorem 12.121 in Section 16.11 

Remark 2.13. The functions Qk and are again Mellin convolutions of qk 
and fk, but now with the Beta density F{1 — t)^~^ supported on [0,1]. We 
will show that P^ in ()2.21l) can alternatively be written as 


Pk{x) 


1 ^ j + ^ + /^)! j 


(2.25) 


ifpk{x) = Thus the factor p,l{n + u)l / {n + u + fi)\ in (|2.24p makes 

Pn{x) monic. The formula (j2.25D is analogous to (|2.16p . and is reminiscent 
of similar formulas in [18], see Remark 12.91 We can also write (|2.23l) as 


Kn(..y) = X-L-‘' 


ds 


27ri \xJ Jy,' s — X Jy t — y 




s — X 


(2.26) 


where S' is a simple counter-clockwise contour encircling both 0 and x. 
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2.4 Generalization 


Proposition Ea Theorem 12.81 Proposition 12.111 and Theorem 12.121 have an 
obvious structural similarity. The new functions in (12.111) and (12.20^ are 
given as Mellin convolutions of fk with a fixed function tp where ip{t) = 
in the case of multiplication with a complex Ginibre matrix as in Proposition 
E71 and p>{t) = F{1 — ^X[o,i](^) ™ the case of multiplication with a 

truncated unitary matrix as in Proposition 12.Ill The two Theorems 12.81 and 
12.121 are obtained as special cases of the following lemma which deals with 
a general function (p. 


Lemma 2.14. Let pQ,... ,pn-i, qo,... ,qn-i be a biorthogonal system for 
the polynomial ensemble D on [0, oo). Suppose 

^A„(7/) det[Ffc_i(yj)]J^^^ (2.27) 

is a polynomial ensemble on [0, oo) with 

Fk{y) = y, y > 0, (2.28) 

where p : [0, oo) ^ [0, oo) is a given non-negative function with finite non¬ 
zero moments. Define 


b^ = 





(2.29) 


(a) Then a biorthogonal system for (|2.27p is given by polynomials Pq, , 
Pn-i, and functions Qo, ..., Qn-i where for k = 0,1,... ,n — 1, 


k k 


Pk{x) = '^aj^kbjX^ if Pk{x) = 

3=0 

^LkX^ 

3=0 

(2.30) 

and 

^ lo (?) 

y>0. 

(2.31) 

(b) Let V’ be given by 



'4){x) = bjX^ 


(2.32) 

where is aiven bv (12.291) for i = 0.n —1 and otherwise arbitrary. 

Assume that the Laurent series (|2.82l) converaes in the annulus \r,i, < 
X < 72^}. Then the polynomial Pk has the alternative representation 


\ ds 
) s ’ 

(2.33) 


where T. is a closed, positive oriented contour lying in the annulus and 
encircling the origin once. 
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(c) Under the same assumptions as in (b), if kn{x,y) is the correlation 
kernel for the polynomial ensemble dni), then the correlation kernel 
Kn{x,y) for the transformed ensemble (|2.27l) is given by 

= (l-f) ■ 

We prove Lemma [2.14l in Section 0] and we obtain Theorems 12.81 and 12.121 
from it as easy consequences. 


2.5 Extensions 

Lemma [2.14l can be applied to other situations as well, in particular to prod¬ 
ucts with a number of Ginibre matrices or a number of truncated unitary 
matrices. We first consider the product of a random matrix X satisfying 
Assumption 11.21 and several complex Ginibre matrices as follows. 

Assumption 2.15. Let ... ,nr be non-negative integers. Let Gi, G 2 , ■ ■ ■, 

Gr be independent Ginibre matrices where Gj is of size (n + Vj) x (n + nj-i) 
for j = l,2,...,r. 

By applying Proposition 12.71 repeatedly, one easily sees that the squared 
singular values of GrGr-i ■ ■ ■ GiX follow a polynomial ensemble, as shown 
in m- It fits in the framework of Lemma 12.141 with the function 


(p = * ■■■ * 


(2.35) 


where (pj{t) = U^e ' for j = 1,..., r, and * denotes the Mellin convolution. 

t ), and we use Meijer 


Then (/? is a Meijer G-function, p{t) = Gq^ 

G-functions freely in our results concerning products with several matrices. 
For their dehnition and properties, see [30], [32|, and [6]. 

Then we have, from Lemma 12.141 


I'l, 


Corollary 2.16. Let Gi,..., Gr be the complex Ginibre matrices from As¬ 
sumption 111.1,71 and let X be the random matrix from Assumption 11.111 with 
£ = n + uq, independent of Gi,... ,Gr- Then the squared singular value 
density of the product Gr ■ ■ ■ GiX is a polynomial ensemble on [0, 00 ) with 
biorthogonal system consisting of polynomials 




-ni,. 


. -nr 


S Pk 


x\ ds 


(2.36) 


where T, is a simple closed counter-clockwise oriented contour around 0, and 
dual functions 


Qkiy') 


^0,r 


ni,, 


t Qk 


y\ (k 
t) t ■ 


(2.37) 
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Furthermore, the correlation kernel for the transformed polynomial ensemble 
is 


Kn{x,y) 


r ° 

2TTi Jy: s Jo t V 0,-j/i,...,-i/,. 


X G 


r,0 

0,r 


ni,...,Ur 




(2.38) 


The proof of Corollary 12.161 will be given in Section [5.21 
Similarly, we consider the product of the random matrix X and several 
truncated unitary matrices as follows. 

Assumption 2.17. Let i'o,i'i,... ,iyr be non-negative integers. Let mj > 
n + nj-i and nij — n — Uj = fij > 1 for all j = 1,2,...,r. Let Ui,... ,Ur be 
independent Haar distributed random unitary matrices of size nij for j = 
1,..., r. Tj are the (n + Vj) x {n + truncations of Uj for j = 1, 2,..., r. 

By applying Proposition l2.11 I reneatedlv. we see that the squared singular 
value density of T^T^-i ■ ■ ■ TiX is a polynomial ensemble, as shown in [25]. 
Then we have the following result as another consequence of Lemma 12.141 


Corollary 2.18. Let Ti,... ,Tr be the truncated unitary matrices from As¬ 
sumption |.2. j7| and X be the random matrix from Assumption 11.111 with 
£ = n + uq, independent of Ti,... ,Tr. Then the squared singular value 
density of the product ■ ■ ■ TiX is a polynomial ensemble on [0, oo) with 
the biorthogonal system consisting of polynomials 


Pk{x) 


1 

2Tii 





0 ,-Vi - /Xi, . . 

0,-lXi,.. 




— Up 



ds 

s 


(2.39) 

where T, is a simple counter-clockwise contour around 0 lying inside the unit 
disk, and dual functions 



Vl+ IJLi,...,Vr+ Lr 


Qk 


y\ M 
t) t ■ 


(2.40) 


Furthermore, the correlation kernel for the transformed polynomial ensemble 
is 


Kn{x,y) 


1 

27rl 




^ '-^r+l,r+l 


0 , fli, . . . , I'r h'r 



QTfi I Vl+ Fl, - ■ ■ + Fr 




Vi,. 


t k„ 


X y 
s ’ t 


(2.41) 


We will prove Corollary 12.181 in Section [6.21 

We remark that Corollaries 12.161 and 12.181 can also be proved by repeated 
application of Theorems 12.81 and 12.121 However, this would lead to cumber¬ 
some manipulations of integrals with Meijer G-functions, which is avoided 
by the use of Lemma 12.141 
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We can use Corollary 12.181 to study the squared singular values of a 
product of the form Ty ■ ■ - Ti. To do so, we can proceed in two ways. We 
can either start with a truncated unitary matrix Ti whose squared singular 
values are a determinantal point process and then apply Corollary 12.181 with 
X = Ti, which is multiplied by r — 1 truncated unitary matrices ■ ■ ■T 2 . In 
this way we recover the result of [251 Corollary 2.6 and Proposition 2.7], as 
we discuss briefly in Remark 16.21 (In Section 15.31 we consider an analogous 
problem, the squared singular values of G,. • • • Gi, in more detail and recover 
the result of [251 Proposition 5.1].) 

The squared singular values of Ti are a determinantal point process only 
if > n, in which case it is the Jacobi ensemble (I2.18P with exponents 
and mi — 2n — ui = fii — n > 0. If /ii < n, then Ti has a singular value 
at 1 of multiplicity > n — ni, and the squared singular values are not a 
determinantal process in the usual, non-degenerate sense. 

However, we may use Corollary 12.181 in a second way, by first taking a 
limit where the random matrix X approaches an n x re deterministic matrix 
with distinct singular values, say 


A = diag(-yau, ^/a 2 , ■ ■ ■, y/a^), CLj G (0, 00 ) are distinct. (2.42) 


Then the squared singular values of T^. ■■ - Ti A are a degenerate form of a 
polynomial ensemble, and we discuss this in Section 15.41 in the context of 
products with complex Ginibre matrices. 

In a further step we can take the limit H —)• /, which in the context 
of Corollary 12.181 leads to the hnal main result of this paper. It answers 
a question posed in |25l Section 7], about the precise conditions on trun¬ 
cated unitary matrices so that the squared singular values of • • • Ti are a 
determinantal point process. 


Theorem 2.19. Let Ti, ..., R. he defined as in Assumption \2. 1 7 with vq = 
0. Then the squared singular values ofT^-'-Ti are a determinantal point 
proeess if and only if 


re < 


(2.43) 


i=i 


If (|2.43li does not hold, then ■ ■ - Ti has a singular value at 1 of mul¬ 
tiplicity > re — Tj ■ 

If (I2.43P holds, then the determinantal point proeess is a polynomial 
ensemble on [0,1] with biorthogonal system consisting of polynomials 


Pk{x) 



U + 1^1 +hi)- j 


A;!G 


0,r-|-l 


A: -I- 1, — rei — Pi,. . . ,-Vr - Pr 
0,-iyi,. . . ,-nr 


(2.44) 
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and dual functions 


Qk{y) = 


J_^r+1,0 

'-^r+l,r+l 



(2.45) 


and the correlation kernel is given by 



pr r(s + 1 + Uj^Vit + 1 + Z^jf + flj) x^y * 

^ r(t +1 + z^)r(s +1 + + /ij) s — t 


.to, —5 — 1 


/^O; • • • ; /^r 

• • • ? 


UX 


(2.46) 


where vq = 0 and ^lq = —n. The contour C is a positively oriented Hankel 
contour in the left-half of the complex s-plane that starts and ends at —oo and 
encircles (—oo,—1], and ^ is a closed contour around [0,re] that is disjoint 
from C. 

We will prove Theorem 12.191 in Section 16.31 

Remark 2.20. In a recent paper [23], Imamura and Sasamoto derived a 
polynomial ensemble related to the O’Connell-Yor directed random poly¬ 
mer model. The biorthogonal system of the polynomial ensemble, see |231 
Formulas (2.1) and (2.3)], has similarities with our formulas (|2.6p and (12.7p . 

Remark 2.21. In this paper we consider the product GrGr-i - ■ ■ GiX or 
TrTr-i ■ ■ ■ TiX with the parameter r hxed. It is possible to consider the 
correlation of the squared singular values of all the products GrGr-i ■ ■ ■ GiX 
or TrTr-i ■ ■ ■ TiX with r = 0,1, 2,... , see the recent paper [33| by Strahov. 


Outline 


In Section (3] we prove Theorem [2T] and Theorem 12.31 about the eigenvalues 
of the sum of a GUE matrix with another Hermitian random matrix. We will 
also discuss some concrete consequences of those results. In Section HI we 
prove Lemma r2.141 In SectionjS] we apply Lemma r2.14l to prove the results in 
Theorem l2.8l on multiplication with a complex Ginibre matrix and corollaries 
on multiplication with several Ginibre matrices. In Section [U] we apply 
Lemma 12.141 to prove the results on multiplication with a truncated unitary 
matrices stated in Theorem 12.121 and the corollaries on multiplication with 
several truncated matrices. At last we prove Theorem 12.191 
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3 Addition of a GUE matrix 

3.1 Proof of Theorem 12.11 


Let be an n X n GUE matrix, and temporarily let M be a fixed n x n 
Hermitian matrix with eigenvalues xi,..., Xn- Then the eigenvalue density 
of H is given by and the random matrix Y = H + M has distribution 




(3.1) 


k=l 


which is the distribution for a random matrix Y in the GUE with external 
source M [9] . It is known from and [21] that the eigenvalues yi,... ,yn 

of y = H + M are distributed according 


Piy,x) = 


1 


(x) 


An(y)det 


Jj,k=l 


y„ = n!(27r)2, (3.2) 


where x = (xi,..., Xn) and y = (yi,..., yn). 

Now, let M be the random matrix satisfying Assumption 1 1.1 1 We obtain 
the joint probability density function of the eigenvalues yi,..., of if + M 
by integrating (13.211 over the eigenvalues xi,..., of M. If their density is 
the polynomial ensemble we obtain that the density of the eigenvalues 

yi,... ,yn oi H + M is (Z„ is defined in (II.ip . the polynomial ensemble of 
the eigenvalues of M) 


P{y) = [ P(y;x)A„(x)det [/fc_i(xj)]" • •-dx, 

JR" 

^n{y) 


ZnZn 


det [/A;_i(xj)]”^^det 




J j,k=l 


dxi 


By the Andreief formula, see e.g. we have 


P{y) = -^^n{y) det 


fk-iix)e-^^y^-P"dx 


dXn. 

(3.3) 

(3.4) 


j,k=l 


with Z'^ = nlZnZm which is the polynomial ensemble defined by (1^ and 

(IT^ . □ 


3.2 Proof of Theorem 12.31 

We assume that pk and for /c = 0,1,..., n — 1, are the biorthogonal system 
as in Assumption 11.11 associated to M. In other words, we assume that for 
/c = 0,..., n — 1, pfc is a polynomial of degree k, is in the linear span of 
fo,...,fn-i, and they satisfy ([13]) • 
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Denote the Weierstrass transform |22] of a function (p by 


1 

y^^{y) = —r^ / dt. (3.5) 

V J—oo 

The integral transform 

1 r+ioo 

/ ^>(s)e2(^-*)'(is (3.6) 

V 27rZ J—ioo 


is the inverse Weierstrass transform: we have (W o = <h for a large 

class of functions. We will need this inversion property only for polyno¬ 
mials <h. A property of the Weierstrass transform is that it sends monic 
polynomials to monic polynomials of the same degree. 

Formulas (j2.6p and (12.7p can now be written as 

Pk{x) = w~^Pk{x), Qkiv) = yVqkiy)- (3.7) 


Since is in the linear span of /o,...,/n-i, Qk is in the linear span of 
yVfo ,..., Wfn-i, and thus also in the linear span of Fq, ..., Fn-i since Fj = 
\f^yVfj by (12.3p . 

By p3.7p . (13.5p and Fubini’s theorem, we obtain 


/ OO /*oo 

Fj{x)Qk{x)dx = / W~^Pj{x)Wqk{x)dx 

-OO J —OO 


/ ^ ^'^^dx'^ qk{t)dt. 


(3.8) 


The expression within parentheses is the Weierstrass transform of W ^Pj 
and is thus equal to pj (t). Hence 


f 


Fj{x)Qk{x)dx = / pj{x)qk{x)dx = Sj^k 


(3.9) 


by (HH). This proves part |(a)| of the theorem. 

Next, we substitute (ESP and (EZP into (ESP and use ESP- This gives 
(|2.8I) and proves part [Tl 

Formula p3.9p is also valid for j = n and fc = 0, l,...,n — 1. Thus Fn is 
a monic polynomial of degree n which is orthogonal to Qoj • • • i Qn-i- This 
implies that Fn is the average characteristic polynomial in the ensemble (12.21) 
and we have proven part (c) of the theorem. □ 
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3.3 Example: GUE plus a unitary invariant random matrix 

In this subsection we assume that the random matrix M from Assumption 
11.11 is defined by the unitary invariant probability measure [TO] 

Cn = (3.10) 

J 

where E is a real-valued function such that the integral defining Cn is con¬ 
vergent and which is allowed to depend on n. This is a random matrix 
ensemble which is invariant under unitary conjugation. If V{x) = x^/2, it is 
simply the GUE. The joint probability density function for the eigenvalues 
of M is the polynomial ensemble with given by fkix) = in (|l.ip . 

For the biorthogonal systems, it is convenient to let pk be the monic degree 
k orthogonal polynomial with respect to the weight on the real line 

and 

Qkix) = ■7-pk{x)e~^^'^\ hk= pk{xfe~^^'^^dx. (3.11) 

n^k J-oo 

The correlation kernel kn can be written as 

n—1 - 

knix,y) = ^ ■j^Pj{x)pj{y), (3.12) 

j=0 

and by the Christoffel-Darboux formula we also have 

^ J_,-V(,)Pn{^)Pn-dy) - Pn{y)Pn-d^) ^ (3 ^3, 

K-i x-y 

Theorem 12.31 applies to this case and we obtain the formula 


Kn{x,y) 


1 Pn{s)Pn-l{t) - pnit)pn-l{s) 

‘27rihyi —J —J — s t 


(3.14) 


for the eigenvalue correlation kernel of the matrix Y = H + M, where H is 
a GUE matrix of size n and M has the distribution (|3.1Up . We show how to 
rewrite this formula in a way that may be suitable for asymptotic analysis. 

For a large class of potentials V (in particular, for n-dependent potentials 
of the form V(x) = nVo(x) with Vq independent of n), large n asymptotics for 
Pn{z) and Pn-i{z) are known for z anywhere in the complex plane. However, 
it is not straightforward to apply saddle point techniques on the integrals in 
(|3.14l) . especially because one has to integrate over the real line, where the 
zeros of the orthogonal polynomials pn, Pn-i are, and where the integrand is 
oscillatory. Therefore, we derive an alternative expression for the correlation 
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kernel, which involves Cauchy transforms of the orthogonal polynomials and 
which avoids integration over the real line. 

Define a 2 x 2 matrix-valued function (which is the solution to the 
Riemann-Hilbert problem for orthogonal polynomials m, IS]) 


Y{z) 


( Pn{z) 

\y-2Trih~^^Pn-iiz) 


^iIZcPnis)e \ 


(3.15) 


The entries in the first column are defined everywhere in the complex plane, 
the Cauchy transforms in the second column are defined for 2 ; G C \ M. 
We can express the correlation kernel kn in (13.1311 in terms of T. It is 
straightforward to check that [lUl Chapter 8] 


i)r-‘Wr(t)Q, (3.16) 

for s G C, t G M. Writing T+(f) (resp. Y-{t)) for the limit of Y{z) as z 
approaches t G M from the upper half plane (resp. lower half plane), we have 
the relation 

y+(t) = y_(t)(^P 1 J’ 

It follows that 

e-'"Wy(t) Q = (y+(t) - y_(t)) (J) , t g m, (3.18) 

and that 

knis^t) = (0 1 ) y-H^) {y+{t) - y-{t)) (?) , (3.19) 

for t G M and for any s G C. In (12.8p . we need to integrate kn{s,t) in t over 
the real line. This integral can be deformed to a contour T which consists 
of a curve in the lower half plane oriented from left to right (for example 
T = M — i6 for some <5 > 0), and its complex conjugate in the upper half 
plane oriented from right to left. We obtain, for any s which is not on T, 


1“ kn{s, t)e-"2(y-^^"dt = ^^ ^ (0 1) y-Hs)y(C) (?) 

if s is inside r, 

+ < (3.20) 

1^0 if s is outside T. 

The last term is the residue contribution from the pole at C = s. In (|2.8|) . 
the integration over the imaginary axis can be changed to any vertical line 
C as contour of integration for s. Let 'w± be the intersection points of C 
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with r, with in the upper half plane and W- = Wj^. (We assume there 
are only two intersection points.) 

Then the contribution from the additional term to Kn{x,y) is 

2m 

= ^ f^{y-x)w+ _ ^{y-x)w-\ ( 3 . 21 ) 

2'iri{ii — x) \ J 

= —^sin {{y - x) Imu;+) , 

7T{y - x) 

where x and y are real. 

We thus obtain the alternative expression for the eigenvalue correlation 
kernel Kn of H + M, 


Kn{x,y) = 


1 


'K{y - x) 


e 2 (^ -y )^{y-x)Kew+ sin{{y - x)luiw+) 


+ 




in which integration over the real line is avoided. This formula may be 
more convenient for asymptotic analysis than (I3.14p . but we do not aim to 
investigate this here. 

Remark 3.1. The random matrix H + M appears in the two-matrix model, 
which is defined as a measure on the space of pairs (Mi, M 2 ) of Hermitian 
n X n matrices, 


(3 23) 


for certain functions Wi and IT 2 such that the above defines a probabil¬ 
ity distribution. If ITi(Mi) = Mf/2, as observed in [13] . the probability 
measure (|3.23D can be written as 


1 e-Tr(i(M3-M2)MlV.(M2)-lM|)^(^^ _ M2)dM2. (3.24) 

This implies that M 2 and Mi — M 2 are independent random matrices: Mi — 
M 2 is a GUE matrix and M 2 is a random matrix from a unitary invariant 
one-matrix model with probability distribution (13.101) . with 

V{M) = W2{M) - ^M^. (3.25) 

In other words, the matrix Mi in the two-matrix model (|3.23p then takes 
the form Mi = H + M, where H is a GUE matrix and M is a random 
matrix from the ensemble (|3.inp . 
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Remark 3.2. If V in ()3.inp is a polynomial, the eigenvalue correlation kernel 
Kn for H + M can be expressed in terms of multiple orthogonal polynomials, 
since we can interpret H + M as, the matrix Mi in the two-matrix model, 
but the asymptotic analysis of those multiple orthogonal polynomials using 
Riemann-Hilbert techniques is hard in general, in particular if the degree 
of V is large. The case where R is a quartic symmetric polynomial has 
been studied in detail in [H] . Our results yield an alternative expression for 
the correlation kernel of the eigenvalues of ff + M, for a general potential 
V, which does not involve multiple orthogonal polynomials but only usual 
orthogonal polynomials and contour integrals thereof. 


4 Proof of Lemma 12.14 


Part (a) From the dehnitions (I2.28p and (I2.3ip it is clear that Qk belongs 
to the linear span of Fq, ... ,Fn-i, since qk belongs to the linear span of 
/o,..., fn-i- It is also clear from (I2.30p that Pk is a polynomial of degree k 
for every fc = 0,..., n — 1. 

To verify the biorthogonality we define an operator £ acting on polyno¬ 
mials p by 


deg(p) 


deg(p) 


Cp{x) = o,jbjX^ where p{x) = 


aiX-' 


(4.1) 


j=0 


j=0 


is arbitrary while bj are defined in (|2.29p . and an operator A4 acting on 
functions q on [0, oo) by 

Mq{y) = (p{t)q y, y > 0. (4.2) 

Then Pk = Cpk and Qk = M.qk by (12.301) and (12.3111 . We prove the identity 

poo poo 

/ Cp{x) Aiq{x) dx = / p{x)q{x) dx (4-3) 

Jo Jo 

and then the biorthogonality follows from 

poo poo 

/ Pj{x)Qk{x) dx = / Pj{x)qk{x)dx = 6j^k (4.4) 

Jo Jo 

since Pj = Cpj, Qk = -Mqk for j, k = 0,1,... ,n — 1. 

To prove (14.3p we calculate by Fubini’s theorem for A: = 0,1,... 

Mqix) dx = <fit) x^q dx^ j. (4.5) 
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We substitute x = tu m. the inner integral and obtain 

poo poo poo 

/ x^M.q{x)dx= / / u^q{u)di 

Jo Jo Jo 


= b 


/ x^q{x)dx. 

Jo 


(4.6) 


By the definition (|4.1D we have Cx^ = b^x^, and thus (14.61) implies that 
(j4.3l) holds if p{x) = for /c = 0,1,... . By linearity (14.311 holds for every 
polynomial p. 


Part (b) The transformation ()2.3n|) is the Hadamard (or termwise) prod¬ 
uct of pk with the function ■i/’- The Hadamard product of two convergent 
Laurent series a{x) = ~ ^ well- 

known (and easy to prove) contour integral representation, namely 


S = «(<.)(.(f)*, 

j=-oa ’’ 


(4.7) 


where T,r is the circle of radius r around the origin with positive orientation. 
Here we assume that the Laurent series for a{x) converges for Va < |x| < Ra, 
and the Laurent series for b{x) converges for < |x| < Ri,. Then formula 
(1321) is valid for ran < k| < RaRb and r G (r^, Ra) n (■^, J^). 

The formula (j2.33p follows from ()4.7p and the definitions (j2.3ni) and 
( 1021 ) . 


Part (c) The correlation kernel for (j2.27p is 


n—1 


Kn{x,y) = '^Pk{x)Qk{y). 


(4.8) 


A:=0 


We insert the integral representations (I2.33P and (j2.31l) for and Qk, and 
interchange the sum with the integrals. Then (|2.34p follows because of 

dop. □ 


5 Multiplication with complex Ginibre matrices 

5.1 Proof of Theorem 12.81 


Theorem 12.81 follows immediately from Lemma 12.141 In the situation of 
Proposition 12.71 we have 

ip{t) = (5.1) 

with v a non-negative integer. Then the moments (I2.29p are 


bj = 


poo 

Uo 


-1 


= [(j + ^^)!] 


-1 


(5.2) 
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and for the function V’ from (j2.32p we take 






(5.3) 


Then parts (a) and |(b)| of Lemma 12.141 give the formulas (12.121) and (j2.13D 
for the biorthogonal system in Theorem 12. lj(a)[ and Theorem E® is also 
derived in the same way. Lemma (2. l^(c)| gives the transformed kernel (12.341) 
in Theorem [22 [b) □ 


5.2 Proof of Corollary 12.161 

Assume X and Gi,..., are as in Corollary 12.161 Applying Proposition 
[27n r times, see also we find that the squared singular value density of 
Y = Gr • • • Gi A is a polynomial ensemble with functions that are the 
Mellin convolution of with 


ip = ipy. * ipr—l * • • • * ip\ 


(5.4) 


where ipj{t) = * for j = 1,... ,r, and * denotes the Mellin convolution. 

Thus is a Meijer G-function 


= gZ 




which has the moments 


nOO J' 

bZ = / t^ip{t)dt = 

k=l 


Now we apply Lemma 12.141 We take for tp 


V’(s) = 'fz 


which is a generalized hyper geometric function 


V’(s) = 


1 


nLi 

and it is also a Meijer G-function 


iPr 


1 

zzi + 1,... , zzj. + 1 


>/>(«) = g!:Li(o_, 


0 


—s 


Then Corollary 12.161 is a direct consequence of Lemma 12.141 


(5.5) 


(5.6) 


(5.7) 


(5.8) 

(5.9) 

□ 
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5.3 Products of Ginibre matrices 

We consider the product 


Yr = GrGr-l ■ ■ ■ Gi 


(5.10) 


of independent complex Ginibre matrices where Gj are defined in Assump¬ 
tion 12.151 with vq = 0. It was shown in [28] that the squared singular value 
density of W is a polynomial ensemble with correlation kernel 


Kn{x,y) 


(27rz)2 jj-^r{u +lyj+ 1) T{v - n + 1) v-u 


(5.11) 

where t'o = 0 and 7 encircles 0 , 1 ,... , re once in positive direction and lies to 
the right of — ^ iM, and T is the Euler Gamma function. For r = 1, this is 
the correlation kernel for the eigenvalues in the complex Wishart Ensemble 
with parameter ui. 

We show how to obtain (15.lip from Theorem 12.81 bv induction on r. 


5.3.1 Base step: Case r = 1 

The r = 1 case of (j5.1ip differs from the well known double contour integral 
formula for the complex Wishart ensemble f |161 Section 5.8]) 


Kn{x,y) 


(2^ Js Jr (u - 1)^ re - u ’ 


(5.12) 


where S and T are disjoint closed contours both oriented counterclockwise 
such that S encloses 0 and T encloses 1. The r = 1 case of (|5.11l) is less 
well known, but it also appeared in the literature as a special case of the 
multiple Laguerre minor process in [H Theorem 3(c)] with the fixed time 
re and special choice of 0:^ = re -|- rei — /c for A; = 1,..., re, see also [ 21 ] for 
its relation to the Laguerre Muttalib-Borodin model. To be precise, the 
relation is 




Knix,y)\^^^ = —K{n,y;n,x) 


(5.13) 


= n + Pi — /c for k = 1 ,..., n 


where K is the correlation kernel defined in [T1 Theorem 3(c)], and the factor 
jg conjugate the kernel into the X]fc=o Pk{x)qk{y) form. The choice 
of contour in [U Theorem 3(c)] is different from that in (|5.1ip . but they are 
equivalent by the residue theorem. In comparison, the better known formula 
(|5.12D is a special case of the Wishart minor process stated in (TJ Theorem 
3(b)], which was originally studied in [ 8 ], [T 2 | and [20|. Both of the two cor¬ 
relation kernels are constructed by Ylk=o Pk{x)qk{y) with pkix) = 
and qk{y) = k\r{k + vi + l)-^ L^I^^\y)y^W~y, by expressing pk{x) and qk{y) 
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into contour integral forms, and then summing them up via a telescoping 
trick. The difference lies in the fact that different contour integral formulas 
for pk and qk are used. 


5.3.2 Induction step 


Let r >2, and assume that (15.lip with r replaced by r — 1 is the correlation 
kernel for the squared singular values of = Gr-i ■ ■ ■ Gi. Then Yr = 
GrYj.-i, and using (j5.11K (with r — 1 instead of r) for kn in the formula 
(j2.14p with = Ur, we obtain a quadruple integral on —1/2 + iM, 7, S, 
and [0, 00 ) for the new correlation kernel. We then first deform the closed 
contour S to a Hankel contour L that comes from —00 in the lower half¬ 
plane, loops around the negative real axis and goes to —00 in the upper 
half-plane. We can then interchange the order of integration. The s- and 
t-integrals are evaluated explicitly as 


1 

27ri 


L 




1 

r(u -|- -i-1) 


= r(u + Ur + 1), 


(5.14) 

(5.15) 


by the integral representations for the Gamma function and its reciprocal. 
Then the result is a double contour integral on — 1/2 -|- iM and 7, that is 
precisely (|5.1ip with parameter r. 


5.4 Limiting case of Corollary 12.161 


We may also approach the calculation of the squared singular values of 
Yr = Gr ■ ■ ■ Gi as a limiting case of Corollary 12.161 

In Corollary 12.161 we assume uq = 0 and let the random matrix X ap¬ 
proach a fixed nx n matrix with distinct squared singular values oi,..., a^, 
or without loss of generality, X ^ A as defined in (12.421) . This is a limiting 
case of a polynomial ensemble with functions fk-i that approach the Dirac 
delta functions 6{x — a^) for every k = 1,... ,n. A limiting biorthogonal 
system is given by the polynomials 

k 

Pkix) = '[\ix - Qj), k = 0,1,... ,n - 1 (5.16) 


with po{x) = 1 and dual functions Qk that are given as a {k + 1) x {k + 1) 
determinant involving Dirac delta functions 


Qkiy) = 


^k+l 


(a) 


1 

Ol 


1 

02 


6{y-ai) 5{y-a2) 


1 

5{y - ttfc+i) 


(5.17) 
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Alternatively, we can use, instead of pk and qk, the Lagrange interpolat¬ 
ing polynomials 

Pk{x) = ^ (5.18) 

and the dual functions %(?/) = 5(y — a^), for A: = 1,..., n. Note that pk has 
degree n — 1 for every k. The biorthogonality 



Pj{x)qk{x)dx = / pj{x)qk{x)dx = 

Jo 


dj,k 


is easy to verify for both systems. 

The correlation kernel for this degenerate ensemble is 


n-l 


kn{x,y) = J 2 Pk{x)qk{y) = '^Pk{x)qj{y) 
k =0 k=l 


n(x-a,) ^ n 


1 


Vi=i 


k=l \j^k 


5{y - cLk) 


Uk — dj j X — dk 


(5.19) 


(5.20) 


and with this expression for kn{x,y) we calculate (|2.38l) by first changing 
variables s ^ x/s, t ^ y/t, and then evaluating the t-integral which results 
in 



The summation can be written as a contour integral, by the residue 
theorem, and we obtain the following result. 


Corollary 5.1. Let Gi,... ,Gr he the Ginibre random matrices defined in 
Assumvtion \2.15\ with vq = 0, and let A be the deterministic matrix defined 
in (12.4211 . Then the correlation kernel of the squared sinqular values of 
GrGr-l ■ ■ ■ GiA is 


Kn{x,y) 




(5.22) 

s — u 
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where T, is a closed, positive oriented contour around 0, and Ca is a closed 
contour, disjoint from S, in the right half-plane encircling each Oj once in 
the positive direction. 

The above derivation of (j5.22p was done under the assumption that the 
Oj’s are mutually distinct. However in (15.2211 we can easily let some or all 
the Oj-’s come together, and the expression is valid for every A with non-zero 
squared singular values ai,..., a„. 

In the limiting case where all Oj —>• 1 and so H = I we obtain 


Kn{x,y) 


1 

{2'KiY 



i 



u 


1,1 

l,r-|-l 


0 




X G 


r,0 

0,r 






{s - 1 

{u — 1)"- s — u 


(5.23) 


with S and T disjoint closed contours where S encloses 0 and T encloses 1. 
This is an alternative expression for (j5.1ip , which for r = 1 reduces to (15.121) . 


6 Multiplication with truncated unitary matrices 

6.1 Proof of Theorem 12.121 


Theorem 12.121 follows from Lemma 12.141 in the same way as Theorem 12.81 
does. 

In the situation of Theorem 12.121 we have 


g}{t) = 


G{1 — for 0 < t < 1, 

0 otherwise, 


with integers u >0 and /r > 1. Then 


6-1 = / ti+^(l - t)>^-^dt = 


-1 . _ (/i- + 


We take for ijj 


= 


1 


]=-V J=-V 


E 


{j + u + p.)\ 

{j + V + ^l)\ 


-X-' 


= pix-^{l - x)-^^-^ 


for 0 < |x| < 1. 


( 6 . 1 ) 


( 6 . 2 ) 


(6.3) 


Then the statements of Theorem 12.121 follow immediately from Lemma 

EH □ 
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6.2 Proof of Corollary 12.181 


The proof is similar to that of Corollary 12.161 Similar to the function in 
(EID, we take here 

(^ = Jf: • • • * 

with now 

^ forO<t<l 

0 otherwise. 


^k{t) = 

Then (p is supported on [0,1] with moments 


=n (7 


(j + i^kV- 


+ 

We can write p as a Meijer G-function 

Vi+ pi,...,l'r + Hr 


p{t) = 

Furthermore, analogous to ()5.7p . we take here 




U + ’^k + f^kV- 


j=0 j=0 \k=l 

which is a generalized hypergeometric series 


(j + t'fc)! 


/ ^ TJ i^k+ fJ-kV- jr 

V’(s) = 11 -—i-r+l-Fr 


fc=l 


l^k'- 


1,1^1 Hi T 1,..., r'j. + /ir + 1 
+ 1 ,..., + 1 


and also a Meijer G-function 

; ! \ /-il^f+l { ~l'l ~ • • • ) ~k'r ~ /^r 

Then Corollary |2.18| is a direct consequence of Lemma [2 


(6.4) 

(6.5) 


( 6 . 6 ) 

(6.7) 

( 6 . 8 ) 

(6.9) 

( 6 . 10 ) 

□ 


6.3 Proof of Theorem 12.191 
6.3.1 Proof of the necessity of (|2.43p 

The necessity of the condition (I2.43p follows from the following linear algebra 
lemma. 


Lemma 6.1. Suppose that 


d = n — Hj ^ 1- 
j=i 


( 6 . 11 ) 


Then the product Yr = Tr- ■ - Ti, where Ti,... ,Tr are defined in Assumption 


2.17 with no = 0, has a singular value at 1 of multiplicity > d. 


27 
























Note that although Ti,...,Tj. are random, Lemma l6.ll uses only their 
matrix structure and the result and argument are deterministic. 

Proof. Each Tj is the truncation of a unitary matrix Uj, which we partition 
as 

=( 1 :) 

where Sj has size fij x (n+Uj-i), and * denotes a block that is not important 
for our present purpose. Thus 


rank(S'j) < Hj. 

(6.13) 

Note that 

(6.14) 

and therefore, since Uj is unitary, 


T*T,+S*Sj = L 

(6.15) 

Below we denote Yj = TjTj^i • • • Ti. Taking j = 1, we 
(16.151) that 

rank(/ — T|*Li) = rank(5('5i) < 

find bv (16.131) and 

(6.16) 


For 1 < j < r we have = Tj^iYj and therefore by (|6.15p 


I - =1- = / - Y*Y, + Y*S*+,S,+,Y,. (6.17) 

Using (|6.13p and elementary properties of the rank, we then get from this 

rank(/ — < rank(/ — Y*Yj) + fij+i. (6.18) 

The two inequalities (j6.16D and (j6.18p . and the assumption (|6.1ip of the 
lemma, lead to 

r 

rank(/ — YfYr) < = n — d (6.19) 

with d > 1. Since YfY^ has size n x n, we obtain 

dimker(/ — YfYr) > d. (6.20) 

Every vector v in the kernel of / — 1^*1^ is a right singular vector of Yj. with 
singular value 1. Thus Yr has at least d singular values at 1. □ 
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6.3.2 Proof of (|2.44p . (I2.45p . and ()2.46l) under condition (|2.43p 

First we consider a limiting degenerate case of Corollary 12.181 with 1^0 = 0 
and the matrix X replaced by the deterministic matrix A defined in (12.421] . 
Then similar to Corollary 15.11 that is proyed in Section [5.41 the biorthogonal 
system and the correlation kernel for the squared singular yalues of T,, • • • Ti^ 
are expressed by (|2.39|) . (|2.4U|) and (|2.41l) with the degenerate form of qk 
and kn containing Dirac 5-functions. The p^s are polynomials and they are 
not degenerate. 

Here we take the biorthogonal system (I5.16p and (|5.17p for the degen¬ 
erate polynomial ensemble of the squared singular yalues of A. Then since 
ip{t)6{y/t — ak)dt/t = ip{y/ak)/ak, it follows from (12.40^ . that 

1 

^/c+1 


p{y/ai)/ai v?(?//a 2 )/a 2 ••• ip{y / at+i) / ak+i 

( 6 . 21 ) 

where p is defined in (16.71) . On the other hand, is obtained simply by 
plugging (I5.16P into (|2.39p . 

Now we take the limit where H —>■ I, so that all aj tend to 1. Then 
becomes the Hadamard product of (|6.inp with 


Qk(jj) 


1 


Afc+i(a) 


1 

ai 


1 

02 


which is 



Pk{x) = Y.i-^P‘ 

j=0 



, _ TT (j +1^1+ fJ-lV- 

' 0- + D)! 


( 6 . 22 ) 


(6.23) 


as also giyen in ()2.44p . 

The limit aj —1 in (|6.21|) has to be done with more care. Note that each 
ipj{x) is zero for x > 1, real analytic on (0,1), with the behayiour Pj{x) ~ 
r(/rj)“^(l—as X ^ 1_. Then inductiyely, for ■ -*^ 1 , 

it is zero for x > 1 and 


p^^\x) = (fij * {pj-l * • • • * <^i)(x) 

(f) *••• *v^i)Wy 


/ 


I' ( 1 (1-.-iK 


c'(l - 


as X ^ 1_, (6.24) 
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for some non-zero constants Cj,c'j. Thus (p = has only + • • • + — 1 

derivatives on (0, oo). 

However, since n < k < n — 1, we have that (p is k times 

differentiable on (0, oo) and then it follows that the limit of (12.401) for all 
ttj —)■ 1 can be evaluated as 


Qk{y) 


(- 

k 



fy\' 

\da y 


a=l 

-p 

a 

Ka) 


(-1)" 

k\ 



(6.25) 


The second equality can be proved by induction on k. Recall that pi is 
the Meijer G-function (16.7D and then by elementary properties of Meijer 
G-functions, we find that Qk is given by (j2.45p . 

Then (I2.44p and (j2.45p give a biorthogonal system for the squared singu¬ 
lar values of Tr • • • Ti, which follow therefore a polynomial ensemble. Having 
the biorthogonal system (|2.44l) and (|2.45p we can proceed as in the proof 
of Proposition 2.7 in |25] to find the expression (12.461) for the correlation 
kernel. □ 

Remark 6.2. In case pi > n, we can also use Theorem 12.121 to prove (j2.46l) 
by induction on r. The argument is similar to the proof of (|5.1ip given in 
Section 15.31 and we do not give the details. We only mention that instead of 
(j5.14l) and (j5.15p we now use 


fJjR 

2'Ki 




i 


1 




r(u ~i~ Pr 
r(/ir)r(u -I- r'r +1) ’ 

r(/ir)r(u + iyr + i) 
Y{v -\- Vj. -\- jx-f -|- 1) 


(6.26) 

(6.27) 


to obtain the desired result. Note that the integral in (I6.26P is equivalent to 
[32l 5.12.9], while (I6.27p is the familiar Beta integral. 

Remark 6.3. For r = 1 the kernel is the correlation kernel of a Jacobi unitary 
ensemble m Section 3.8.3], [M!- The double contour integral formula (|2.46l) 
with r = 1 for this kernel is not well-known, and it was first described, to 
the best knowledge of the authors, in [U Theorem 3(d)] as a special case 
of the Jacobi-Pineiro minor process. (See also [2T] for its relation to the 
Jacobi Muttalib-Borodin model.) 

In the notation of [T], we take fixed time n, = n + ui — k, for k = 
1,..., n and M' = pi. Then 


Kn{x,y)\^^^ = 


/^1 


-K{n,y,n,x) 


= yi ui — k for k = 1,... ,n, M' = fii 


(6.28) 


analogous to (|5.13l) , where K is the correlation kernel defined in [TJ Theorem 
3(d)]. 
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